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ABSTRACT 



A layer-ocean model was developed based on the mixed 
layer models of Kraus and Turner (1967) and Denman (1972), 
The ocean model was coupled to Pearson’s (1972) time vary* 
ing, symmetrical-stationary hurricane model which was 
based on a model proposed by Riehl (1963), and in which 
the air-sea interaction was specified by Cardone 1 s (1969) 
extension of Blackadar's (1965) two-layer, baroclinic 
boundary- lay er model. Time-dependent solutions for ocean 
mixed layer depth and temperature were obtained in res- 
ponse to interaction with the atmospheric model. Solu- 
tions indicated that the interaction between entrainment 
mixing and upwelling was most important in changing mixed- 
layer depth and temperature. Radiational effects, inter- 
nal waves, turbulent scale energy dissipation, and large- 
scale ocean currents were not included in the present 
ocean model . 
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INTRODUCTION 



I . 

It is widely accepted that the source region for the 

hurricane heat engine is the underlying ocean. Although 

the ocean is virtually an infinite heat source, it is 

generally accepted that sea-surface temperatures in excess 

of 26C are necessary for forming and maintaining a mature 

hurricane. At the same time the ocean acts as a sink for 

the hurricane's momentum which is imparted to the ocean 

through surface stress. It is the balance of these, two 

processes which was the basis for Riehl's (1963) steady- 

state hurricane model in which the tangential wind profile, 

outside the radius of maximum winds, was specified as 

1/2 

Vgr = constant. 

Corgnati (1971) adopted Riehl's work to develop a 
model which used boundary- 1 ay er predictive equations for 
moisture and temperature, and bulk aerodynamic equations 
for the fluxes at the air-sea interface. Corgnati in- 
cluded a fairly simple ocean model which described the 
oceanic response to convective mixing as a result of heat 
extraction; however, most experiments were run with con- 
stant sea surface temperatures. 

Pearson (1972) modified Corgnati's model by using 
Cardone's (1969) extension of Blackadar's (1965) two- 
layer near- neu tra 1 , baroclinic boundary -1 ay er model to 
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calculate the fluxes at the air-sea interface. Pearson's 
model produced steady-state solutions for constant sea- 
surface temperatures, and exhibited greater sensitivity 
to these temperatures than did the Corgnati model. 

It was felt that to more realistically describe the 
hurricane-ocean interaction, a more complex ocean model 
incorporating wind mixing, convective mixing, and advec- 
tion, was required. It was the object of this study to 
develop such a model and couple it to Pearson's atmos- 
pheric model through the boundary-layer fluxes of heat 
and momentum. The current ocean model is based on mixed- 
layer models of Kraus and Turner (1967) , which predicts 
the mixed-layer temperature and depth in response to wind 
mixing and convection, and Denman (1972), which includes 
the effect of upwelling. It was hoped that the ocean 
model may help to explain which of the mechanisms is of 
greatest importance in producing the areas of low sea- 
surface temperatures observed in the wake of hurricanes. 
Investigators have attributed the low temperatures to 
various mechanisms with little agreement as to which is 
the most important. For example, Jordan (1964) concluded 
that mixing was the more important factor. Fischer (1958) 
considered upwelling to be the most important. Leipper 
(1967) concluded that upwelling, produced by surface 
divergence, was the most important near the hurricane's 
center and that mechanical mixing dominated at larger 
rad i i . 
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The present ocean model was run using each of the 
three cooling mechanisms (convection, wind mixing, and 
advection), separately and in various combinations, in 
an attempt to study their relative effects upon the 
ocean upper thermal structure and subsequent heat content. 
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II. DESCRIPTION OF THE MODEL 



A. ATMOSPHERIC VORTEX 

Riehl (1963), in his early work, assumed conservation 
of potential vorticity in the inflow layer and a bulk 
aerodynamic equation for surface stress with constant drag 
coefficient and inflow angle, to derive the relation 



i / 2 

v Q r = constant 

W 



( 1 ) 



which specifies the radial profile of cyclonic tangential 
winds at the top of the atmospheric boundary layer as 
shown in Figure (1). The outer limit of the circulation 
described by Riehl 1 s model was defined by the radius (r^) 
at which the cyclonic flow becomes anticyclonic in the 
outflow layer. Riehl assumed that the air ascending in 
the eye-wall region, and flowing cyclonically outward 
aloft, conserves absolute angular momentum. With this 
assumption, one obtains 



r 

o 







1/2 



(2) 



where f is the Coriolis parameter and the i subscripts 

denote values at the radius of maximum wind speed. The 

present atmospheric model assumed that the radius of 

maximum wind (r.) was coincident with the eye-wall radius. 

i 

Since the eye wall is a preferred region for deep convec- 
tion and latent heat release in observed hurricanes, 
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those boundary- lay er air parcels with the greatest heat 



and moisture content in the model were assumed to be at 
the eye wall. Thus, r_^ was taken to be the radius at 
which the equivalent potential temperature (0^), averaged 
over five grid increments (Ar = 3km) , was a maximum. 



Between the vortex center and r. the winds were assumed 

1 

to obey v r ^ = constant, which agreed with data from 
o 

Riehl (1963), as well as Shea and Gray (1972). 

From observations, Riehl developed the relationship 



P ' = -2 . 56 0 1 

e 



(3) 



where P 1 is the pressure departure from 1005 mb and 0 ' 

is the equivalent potential temperature departure from 
350 K. It was assumed that the atmospheric vortex was 
in a state of qua s i - equ i 1 ibr ium with the boundary layer 
fluxes of heat and moisture, which determined the in- 
crease in 0^' and thereby P* of the inflowing air. 

Assuming a gradient balance between the wind and 
pressure field 



I UL 

p 3 r 




f v 9 



(4) 



and using the relation 



v . 




(5) 
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integration of (4) from r^ to r^ , yields the following 
expression for the maximum wind speed: 

0 . 5r 



v 



0 . 

i 




[ - 2 f /r . r 
1 o 



r . 
i 



+ /4f 2 (/r . r - 
1 o 



r . ) 2 + 10 6 r -r . (0 

i o i e . 

l 




( 6 ) 



For the symmetric model, the time-dependent predictive 
equations for potential temperature (0) and specific 
humidity (q), which characterize the boundary- 1 ay er pro- 
perties through a specified depth (AP), are 



9j0 

at 



3 (rv 0 ) 
r 9 r 



a ( 030 ) 

ap 



gQ, 



c Ap 
p 



and 



~ a (rv q) ^ , x gQ 

9q r _ 3 (cog ) 

at r ar ap lap 



(7) 



( 8 ) 



where (Q s ) an d (Q E ) are the oceanic sensible and latent 
heat fluxes, respectively. and L are the specific heat 

and latent heat of vaporization. Radial (v ) a nd vertical 
(U)) components of motion are related through the con- 
tinuity equation 



a co 

ap 



9(rv r ) 
r 9 r ' 



(9) 



The surface fluxes Q , Q , T n , and the integrated radial 

s e e s 

motion ( v r ) are determined from the boundary-layer model. 
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B. BOUNDARY-LAYER MODEL 



Following Pearson (1972), the boundary- lay er model 
was a two-layer, baroclinic model developed by Cardone 
(1969) as shown in Figure (2). The Cardone model of the 
marine boundary layer was an extension of the near- 
neutral, f ixed- ter rain model described by Blackadar (1965). 

The model used the v Q profile from the atmospheric vortex 

u 

model, together with a computed boundary- lay e r thermal 
wind and made use of Monin-Obuhkov similarity theory and 
stationary Ekman-layer theory to calculate the friction 
velocity u*, the stability index L*, and the inflow angle 
a. These are then used to calculate a vertical mean radial 
wind, u, within the surface layer from 



u 




[Ln ^ - ip (L*) ] J / 4 



z . 



o 






( 10 ) 



and the total heat flux (H) from 



H = Q s + Q e = [ -u *C^pa^ 0 ] /KgL ' 



( 11 ) 



where 



1/2 

u*= friction velocity, (T/p) , where X is 

the radial stress component. 



Z 



o 



6.84 x 10 5 +4.28 x 10 3 u* -4.43 x 10 



surface roughness (m) , 



i|;(L^)= wind profile function [Cardone (1969)] 



L = Z/L 1 , stability index 



Z_^= heights in the surface layer 
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K = von Karman constant, 

0 = mean potential temperature in the 
boundary layer, and 

a^= K^/K^, ratio of heat transfer coeffi- 
cient to turbulent transfer 
coefficient for momentum. 

The individual amd fluxes were specified accord- 

ing to the Bowen ratio C (0 -Q ) /L (q -q ) . The radial 
^ p a x w v ^a ^w 

wind within the spiral layer was determined through in- 
tegration of the Ekman solution to get 



/ V dZ = /2 V sin a 

, r 0 2Bh 

n 



-Bh 3-tt 

[ cos (a+ — -Bh) 



- sin (a+ 



3tt 



“ Bh ) ] 



( 12 ) 



where 



B = ( f / 2u * Kh ) 



1/2 



To include the* centripetal acceleration, the Coriolis 

v © 

parameter (f) was replaced by (f + ) . The values of 

v^ in the surface layer (Equation 10) and Ekman layer 
(Equation 12) were then used in determining the vertical 
motion (Equation 9) and the advection of heat and 
moisture [Equations (7) and (8)] in the prediction 
portion of the model. 



C. OCEAN MODEL 

1 . De s cr ipt i on 

In an attempt to more realistically describe the 
air-sea interaction associated with a time-varying 
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hurricane model, it was decided to develop an ocean model 
capable of producing time-varying solutions of mixed- 
layer depth and temperature in response to atmospheric 
forcing. Previously the hurricane model produced time- 
varying solutions for constant sea-surface temperatures 
until a steady state was reached. 

The present ocean model was initialized with sea- 
surface temperature (T^) / mixed-layer depth (h^) / below 
layer gradient (3T/3Z), and deep-layer temperature (T^) 
as shown in Figure (3). By allowing the ocean thermal 
structure to change, it was hoped that a more thorough 
understanding of the air-sea interaction process might 
be possible. The present model is designed for a wind- 
dominated region where there is negligible downward heat 
flux, either by solar radiation or sensible and latent 
heat (condensation) fluxes from the atmosphere to the 
ocean. In the absence of upwelling this precludes any 
possibility for a decreasing mixed-layer depth. The 
atmospheric forcing of the ocean model was by surface 
kinetic energy input (wind stress) and upward sensible 

and latent heat (evaporation) fluxes. Both wind stress 

/ 

and heat flux produced changes in mixed-layer depth and 
temperature, surface stress by entrainment mixing and the 
heat fluxes by convective mixing. The entrainment mixing 
was caused by surface generated turbulent motions at the 
depth of the mixed-layer interface. 
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A partitioning procedure w as used to budget the 
total atmospheric stress (T^) between surface turbulence 
production and generation of a radial mass transport ac- 
cording to Ekman theory. The effects of entrainment 
mixing (turbulence) were assumed to be redistributed 
uniformly throughout the layer by turbulent diffusion. 

A radial wind-driven current was computed from the mass 
transport and, by employing the continuity equation, areas 
of upwelling and downwelling were calculated. Thus, the 
ocean model incorporated the primary mechanisms believed 
to produce the observed ocean cooling in the wake of 
hurricanes. Internal waves, large-scale currents, and 
turbulent-scale energy dissipation within the mixed layer 
were not included. At six-hour intervals the ocean heat 
content was computed and compared to the initial values 
as a check on the model's solution of mixed-layer depth 
and temperature. 

The time-dependent mixed-layer temperatures 
generated by the ocean model produced variations in the 
atmospheric model's v Q profile, which in turn produced 
variations in the atmospheric forcing (surface stress and 
heat fluxes). This tended to bring the coupled models 
into a quasi-steady state condition in all cases except 
ones in which advection was included. 

2 . Mathematical Development of Ocean Model 

By neglecting molecular heat fluxes and viscous 
generation of heat, Denman (1972) wrote the first law 
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of thermodynamics (conservation of thermal energy) for 
an incompressible fluid in the following form: 



dT 

dt 



q t 

P c 

w p 



(13) 



where 

= y R * e 1 #■ specifies the heat source terra, (14) 

p = ocean water density, 
w 

y = average extinction coefficient, 

0^ = specific heat at constant pressure, 

R* = solar radiation incident on the sea surface, 
and 

Z = depth in the ocean, positive downward. 
Substitution of (14) into (13) gives the first law in 
the form 



dT 

dt 



YRe 



-yz 



(15) 



where 



R = R*/PC 
* P 



The total derivative was expanded in flux form and (15) 
was written in time - aver aged turbulent form: 



or 



3T 1 9 (u r rT) 9 (w T) 3 (w ' T ' ) 
3 1 r 3r 3 z + 3z 



YRe 



-YZ 



(16) 



dT _ -YZ 3 (w- T' ) 

dt " YRe 3z 



(17) 



where u^ = vertical average radial velocity in the mixed 

layer, and T, u and w were replaced by T + T 1 , u + u 1 , 

r r r 
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and w + w', respectively. The term involving (u^T 1 ) in 
Equation (16) was dropped because it was assumed that 
there was no systematic correlation between u* and T*. 
Equation (17) is the heat conservation equation where 
9 (w , T , )/ 3Z represents the local divergence of the turbu- 
lent vertical heat flux. It was this term which redis- 
tributed the heat exchanged at the boundaries uniformly 
throughout the homogeneous mixed layer. At the ocean 

surface, the downward turbulent heat flux, (w'T') , was 

o 

equal to the nonradiative heat transfer through the 
ocean surface, i.e., 



( w 1 T 1 ) 

o 



= - (Q + Q ) /p C 

^E *S 7 H w p 



(18) 



where Q and Q were the upward latent and sensible heat 

b b 

fluxes, respectively, computed in the boundary -layer 
model. At the bottom of the mixed layer, the downward 
turbulent heat flux was equal to the heat transfer through 
the interface due to entrainment mixing, i.e.. 



( w 1 T ' ) h 

o 



dh 




W) { V T h } 



(19) 



where dh 

o 

— — = time variation of mixed-layer depth, 
d t 

w = vertical advection velocity at depth (h^) / 

T = mixed-layer temperature, and 

= temperature at a finite depth below mixed- 
layer depth. 



For entrainment mixing, the turbulent flux (w'T 1 ), was 

h 

o 

positive which corresponds to warm water (T*> 0) being 
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forced downward (w*> 0) and cold water (T'< 0) being 
lifted upward (w'< 0). This process was accomplished 
against buoyancy forces by the kinetic energy input at 
the surface and the potential to kinetic energy conver- 
sion by convective overturning. 

Integrating Equation (17) over the depth of the 
mixed layer (h^) and replacing the boundary turbulent 
heat fluxes by Equations (18) and (19), yields a thermal 
energy equation for the mixed layer, 

dT 

ar h o - R(1_e " Th) " is e +2 s )/p w c p ~ ( ftT ( VV <20 > 



By expanding the total derivatives in flux form and ig- 
noring the radiation effects, (R), (20) becomes a predic- 

tive equation for the local change in the mixed-layer 

depth (h ) , i . e . , 
o 



9h 
• o 

9t 




3T s 1 3(u r rT s ) 3 (w T s ) 
3t + r 3r + 3Z 



^E + 

p w c p ( VV 



- u 



3h 

c 

77 



4- w 



( 21 ) 



To develop a corresponding equation for mixed-layer tem- 
perature changes, Denman (1972) integrated Equation (17) 
twice, once over a depth z and once over the depth of the 
mixed layer (h^) . The resulting expression, 
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dZ 



( 22 ) 



dT 



S 



dt 




( v y 

P „ C p 



h 

o 



h 

o 

/ ( w ' T 1 ) 

Z 

o 



represents a mechanical energy equation for the mixed 

h 

layer. The term / (w'T')dZ specifies the conversion 

o 

of potential energy into kinetic energy by convection 
within the mixed layer. If one considers the mechanical 
energy balance expression 



W + G - D = 0 



(23) 



where G = kinetic input from the wind, 

D = dissipation within the mixed layer, and 
h 

W ="/ ° ( w * T *) dZ (potent ia 1 to kinetic energy 

o conversion ) , 



substitution of Equation (22) into (23) yields 



dT h 2 
S o 

dt 2 



<V V 

P w c p 



h 

o 



(G-D) . 



(24) 



Expanding the total derivative of temperature in flux 
form and solving for the local temperature change yields 




9t 



2 (G-D) 



2 <Q E - 

h o P 



Q s } 

C 

w p 



1 

r 



3 (u r rT s ) 
3r 



8 (w T g ) 
3 Z 



(25) 



which is a mechanical energy equation specifying the 
change in mixed-layer temperature due to entrainment 
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mixing, convection, and advection. Equation (25) can 

be used to eliminate 3 t /3t from the heat equation (21) : 

S 



9h o _ 2 (G-D) 

3t h (T -T, ) 
o s h 



(q e* V 

P C (T -T ) 

w p s n 



u 

r 



3h 

o 

3r 



+ 



w 



(26) 



Equations (25) and (26) represent a system of two equa- 
tions with three unknowns, T , h , and T, (u and w will 

so hr 

be obtained from Ekman theory and the equation of con- 
tinuity). To close this system therefore, a third equa- 
tion or relation was required which specified changes in 
T^ as being dependent upon changes in h . The parameter 
which imposed this dependency was the below-layer gradient 
(3T/3Z) which was assumed to be uniformly affected by ver- 
tical velocities and thus remained constant. Time varia- 
tions of T were specified by the expression 



flh 

3t 



- (w - 



dh 

~) ( — ) 

dt ; K dt } h 



(27) 



where w = vertical advection velocity. 

The solution for T^ , which represented the lower tempera- 
ture of the entrained water, proved to be very crucial 
in the mixed-layer depth and temperature solutions. It 

was the temperature difference (T - T ) which imposed the 

s n 

primary retarding force against unlimited deepening. As 
time progressed the difference (T^- T^) increased, forming 
an actual step in the temperature profile at the base of 
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the mixed layer. Because of the stable lapse-rate 

(8T/9Z < 0), upwelling (w<0) contributed to a lowering 

of T, and, hence, increased the difference (T - T, ) . 
h s h 

Thus, upwelling had the effect of retarding the tendency 
for mixed-layer deepening. At the same time upwelling 
(w<0) enhanced mixed- layer temperature changes by de- 
creasing the layer thickness (Equation 26) so that the 
energy input and heat fluxes would have a greater effect 
(Equation 25). Thus, Equations (25) (26), and (27) 

formed a system for unknowns h Q , T^ , and T^ . Space and 
time derivatives in (25) and (26) were approximated by 
centered differences, except for the initial time step 
which used a forward difference. A forward step was 
also used every 17th time step (3 hours) to prevent di- 
verging solutions. Forward time differences were used 
in (27). The difference forms of (26) and (27) were 
solved together for h^ and T^ using a first guess depth 
and temperature (T^) from the previous time step. The 
first guess depth in (26) and (27) was increased by 
0.05m increments until the difference between the right 
and left hand sides of the difference form of (26) was 
less than or equal to 0.1m. Using the new mixed-layer 
depth, mixed-layer temperature was calculated from (25). 

The mechanical energy input (G) , in the ocean 
predictive equations (25) and (26), was specified 
according to Kraus and Turner (1965) as 
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G 



1 



(28) 



gotp . 



/a, 3/2 
(— ) u t 



w 



where g = gravitational acceleration, 

1 dp w 

a = — — , coefficient of expansion, 

P w 

p = air density, 
a 

p = sea water density, 
w 

u* = atmospheric friction velocity, and 
= atmospheric surface stress. 

In the ocean model, surface stress (T ) used for the 

o 

energy input term (G) depended exponentially upon the 
ratio of the mixed-layer depth (h^) and the roughness 
leng th ( ) : 



T 

o 



T EXP 
a 




(29) 



where 



T = total stress available from boundary “lay er 
model , and 

C = an arbitrary constant. 



h 

For a given roughness length (Z ) , the function exp (- *r — ) 

O C Z 

o 

specified the percentage of total stress (T ) ( available 

cl 

for turbulence generation as a function of mixed -layer 
depth (Table I). The use of the above partitioning pro- 
cedure was an attempt to budget the available energy 
between turbulence generation (G) and production of wind- 
driven currents in the ocean model. Experiments were run 
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with various values of C and compared to each other as 

well as to a run in which all the surface stress was 

utilized for turbulence generation. The dissipation, D, 

was set to zero in all experiments. 

The radial velocity, u^, which is needed to close 

the system, was calculated in accordance with Ekman theory 

which specifies a spiral layer to a depth (Z ) where the 

e 

current direction has reversed from its surface direction 
and the speed is exp (-7T) of the surface speed. Below 
the Ekman layer it was assumed that the wind had no direct 
influence. According to Ekman the radial component of 
motion within the spiral layer was described by 



(T -t ) cos a 

u = T-y EXP (-BZ) [cos (-BZ) - sin(-BZ)] (30) 

p ( 2 f ) 7 

w p 

1 t-T 



where 



Z^ = 7T/B, specifies the Ekman depth. 



B = [ 



f j 1/2 



2y/p 



w 



f = coriolis parameter, 

= sea water density, 

y = vertical eddy coefficient of viscosity, 

Z = depth, and 

a = atmospheric inflow angle. 

Integration of (30) over a specified Ekman depth (Z^) 
yields an expression for the net mass transport in the 
radial direction for the Ekman layer, 



M r = V f 



(31) 
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TABLE I 



Stress used for mixing (T ) vs. 
mixed-layer depth as a % of 

total stress (T ) . Z =10 2 m, 

C =10 and 5 x a l0 . ° 



Mixed layer 
Depth (m) 

5 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 



T as % of x 
o , _ , ^4 x a 



.95 
. 90 
.82 
. 74 
.67 
. 61 
. 55 
. 50 
.45 
.41 
.37 



T as % o f x 
o . _ _ _ _ 3 , a 

(C = 5x 10 v ) 

. 90 
.82 
. 67 
. 55 
.45 
.37 
.30 
.25 
. 20 
. 17 
. 14 
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where 



= (l - T ) cos a, tangential atmospheric surface 
U a o 

stress. 

This represents a radial mass transport per unit distance 
in the tangential direction. One can calculate a mean 
radial advection velocity (u^) f° r the mixed layer from 
the relation 



u Z = M /p . (32) 

re r w 

The mean radial (u^_) and vertical (w) components of motion 
are related through the continuity equation 



9 ( P W w) = _l 9 ( P W u r r) 

3 Z r 3 r 



(33) 



Equation (33) was integrated over depth (Z^) using the 
boundary condition w(z=0) = 0, to obtain w at Z 



i3 , r T e 

P w W = - r 37 (_ f 



(34) 



which was taken to be w at h . The probable occurrence of 
a return flow at some depth below the mixed layer was not 
included in the present ocean model. Thus, the only effect 
on the below“layer gradient (3T/3Z) was that due to verti- 
cal velocities and, hence, the gradient remained constant. 

As a check on the cooling processes (wind mixing, 

convective mixing, and advection), the initial ocean heat 

2 

content (cals/cm ) for a depth Z, = 130m was calculated 

D 

based on the initial temperature profile (Figure 3). This 
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value was then compared to the heat content at each grid 
point, computed at 6-hour intervals. For convective mix- 
ing-only, with the kinetic energy input (G) and the 
advection velocities (u^ and w) equal to zero, Equations 
(25) and (26) predicted excessive mixed-layer deepening 
and cooling. Also the final grid point heat contents did 
not reflect just the heat (H) lost to the atmosphere due 

to upward heat fluxes (Q and Q ) . Hence, for this case 

E S 

the mixed-layer depth and temperature were specified by a 
thermal equilibrium problem using the following set of 
equations : 



where 



T -T 1 

h T = h ' (T QAt/p C h')+ Ah • % S 
os o s wpo 2 



h 

o 

T 

s 



h 1 

o 

T 1 

s 



+ Ah 




(35) 



h ' = 
o 

T 1 = 
s 

8T/3Z. = 
Q = 




initial mixed-layer depth, 

initial mixed-layer temperature, 

be low- layer gradient , 

Q + Q , total upward heat flux, 

E S 

sea water density, 
time step (sec) , and 
specific heat 



For convective mixing only, equations (35) were solved at 
each time step for the change in mixed-layer depth (Ah) . 
With this calculated depth change, the new layer depth (h^) 
and temperature (T^) were specified. 
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Figure 2. Boundary layer model. 
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Figure 3. Ocean model. 
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III. DISCUSSION OF RESULTS 



As stated previously, the object of this study was to 
develop an ocean model, capable of producing time-dependent 
solutions for mixed-layer depth and temperature in response 
to forcing from a symmetrical, stationary hurricane model. 
Previous work with the hurricane model was done with con- 
stant sea-surface temperatures and, consequently, little 
was learned of the ocean thermal changes or of the corres- 
ponding changes produced in the hurricane model. 

In this investigation three mechanisms believed to con- 
tribute most in changing the ocean thermal structure were 
isolated and studied. These were entrainment mixing across 
a stable layer (wind stress), convective mixing (latent 
and sensible heat fluxes), and vertical and horizontal 
advection. It was hoped that the combined time-dependent, 
atmosphere-ocean model would help determine which of these 
mechanisms, either acting alone or in combination, was most 
important in producing the cold sea-surface temperatures 
observed beneath hurricanes. To accomplish this, a series 
of experiments were run using each of the cooling mechanisms 
alone and in combination with each other. Figures (4) 
through (32) depict radial profiles of various ocean para- 
meters for the different experiments. For these runs the 
ocean model was initialized with mixed-layer depth of 
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30 meters, mixed-layer temperature of 30C, below-layer 

gradient of lOC/lOOm, and deep-layer temperature of 20C. 

The lower boundary of the model was fixed at 130 meters 

(Figure 3) f below which there were assumed to be no 

thermal effects. Other initial temperature profiles can 

be used; however, the above values were believed to be 

fairly representative of the tropical ocean regions. The 

atmospheric model was initialized with air temperature of 

29C, radius of maximum wind (r_^ = 27km), and maximum wind 

(v A = 30m/sec). In runs with various constant sea-surface 
0 . 
i 

temperatures, Pearson (1972) found a steady-state Vg of 

i 

60m/sec for 30C. 

The value of the constant (C) in the wind stress parti- 

4 

tioning equation (29) was arbitrary fixed at 10 in all 

the initial runs except the advecti on-only case. For this 

run the partitioning process was not used; consequently, 

all the available stress (T ) was used to produce radial 

a 

wind-driven currents (u ) and vertical velocities (w) . 

r 

Since the partitioning process was an important assumption 
in the model, an entrainment-convec t ion-advect ion experi- 
ment was run with a value of 5 x 10^ for the constant (C) 
in Equation (29). With a smaller value of C, Equation (29) 
partitioned more atmospheric stress to producing ocean 
currents and less to entrainment mixing at all mixed-layer 
depths (Table I). An entrainment-convec t'ion run was also 
made in which the partitioning process was not used, hence 
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all available stress from the boundary-layer model was 
used in the kinetic energy term (G) in the mixed-layer 
depth and temperature predictive Equations [(25) and (26)]. 

A. CONVECTIVE MIXING ONLY 

For the convective mixing-only case, Equations (35) 
were solved for the mixed-layer depth (Figure 4) and tem- 
perature (Figure 5). The atmospheric forcing was by 
latent and sensible heat fluxes from the ocean surface to 
the atmosphere as computed in the boundary -lay er model. 
Maximum cooling (0.5C) and deepening (5 meters) of the 
mixed layer at 18 hours occurred under the region of maxi- 
mum wind (r^) as expected from Equation (11) which relates 
total oceanic heat loss (H) to the cube of the friction 
velocity (u*). Oceanic heat loss computed from initial 
and final temperature-depth profiles (Figure 6) in the 
ocean model compared quite favorably with the total heat 
extraction (H) computed in the boundary- lay er model 
(Figure 7). The amount of ocean cooling due to convection 
was comparable to Jordan's (1964) estimates. 

B. ENTRAINMENT MIXING 

Radial profiles of mixed-layer depth and temperature 
due only to entrainment mixing (Figures 8 and 9) were ob- 
tained by iterative solutions of Equations (25), (26), 

and (27). There was no ocean current generation in this 

4 

run, and the partitioned stress (29), with C = 10 , was 
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Figure 4. Radial profile of mixed-layer depth at 
18 hours due to convective mixing. 




Figure 5. Same as Figure 4 except mixed-layer 
temperature . 
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Figure 6. Predicted 18-hour vertical temperature 
profile at 30 kmfor convective mixing. 
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Figure 7. Predicted 18-hour radial profiles of 

ocean heat loss and accumulated (Q E +Q S ) 
from boundary layer model. 
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Figure 8. Radial profiles of mixed-layer depth at 
6 and 18 hours for entrainment mixing. 




Figure 9. Same as Figure 8 except mixed-layer 
temperature . 
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the forcing which produced entrainment mixing through sur- 
face turbulent generation (G) . Maximum cooling (0.5C) 
and deepening (23 meters) at six hours and 1 . 5C and 45 
meters at 18 hours again occurred beneath the region of 
maximum kinetic energy input (28) as expected from (25) 
and (26). The increased cooling and deepening w ith time 
at larger radial distances was due to the outward migra- 
tion of the hurricane eye wall (r^) in response to colder 
sea-surface temperatures near the vortex center. 

Since entrainment mixing produced only vertical heat 
redistribution, the heat content of the ocean beneath the 
atmospheric vortex should be conserved except for numeri- 
cal errors. Figure (10) depicts the radial distribution 
of oceanic heat loss at 18 hours for the entrainment-only 

experiment. The area averaged heat loss for the solution 

2 

domain (Table II) was 137 cals/cm . Considering the 

sensitivity of heat content to changes in temperature and 

depth and since the depth and temperature solutions were 

computed numerically, and are subject to truncation errors, 

the values in Figure (10) are considered within acceptable 

limits. For example, in the region of maximum mixed-layer 

depth (75m), a temperature difference of 0.1C will give 

2 

a heat contant change of 750 cals/cm . The maximum error 

2 

in Figure (10) is 800 cals/cm . 

Figure ,(11) depicts ocean temperature profiles at 6, 

12, and 18 hours as computed in the model. Note the strong 
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1000 T 




Figure 10. Predicted 18-hour radial profile of 

ocean heat loss for entrainment mixing. 
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Figure 11. Predicted 6- , 12-, and 18-hour vertical 

temperature profiles at 30 ^m for 
entrainment mixing. 
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gradient which develops below the mixed layer. As men- 
tioned in a previous section, the temperature difference 

(T^- T ) is the primary retarding force which balances 

3 T 

the atmospheric forcing. The below-layer gradient (7: — ) 

o Z 

is an important factor in determining the depth at which 
this balance will occur. For a given energy input (G), 
one expects deeper mixing for a weak gradient and 
shallower mixing for a strong gradient. 

C. ADVECTION-ONLY CASE 

For the advec tion-only case (G=0, and Q + Q =0), mixed” 

E S 

layer depth and temperature changes were determined by the 

advection terms in the predictive equations [(25) and (26)] 

The partitioning process (29) was not used, hence, all 

available atmospheric stress (T ) was used to produce a 

a 

radial mass transport (32) and a mean radial wind-driven 
current (33). -The rapid increase in u^ (Figure 12) from 
the vortex center to the radius of maximum wind (r^) pro- 
duced a region of horizontal divergence and, hence, enor- 
mously strong upwelling (Figure 13). Just beyond radius 
a small region of convergence and, hence, downwelling 
occurred. However, the gradual decrease in u^ was 
balanced by the increased radius to maintain a fairly con- 
stant total mass transport. Consequently, areas of signi- 
ficant downwelling were not evident within the solution 
domain (300* km) . 

The effects of strong upwelling were quite evident in 
the mixed-layer depth and temperature profiles (Figures 
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Figure 12. 6 - and 18-hour radial profiles of ocean 

radial velocity for advec tion-only (no 
partitioning ) . 




Figure 13. Same as Figure 12 except ocean vertical 
velocity . 
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14 and 15). Within six hours the mixed layer, near the 
vortex center, was eliminated and the sea-surface tempera- 
ture decreased until the cold deep ocean water (20C) was 
upwelled to the surface. At 18 hours the radius of maxi- 
mum wind (r^) had moved outward in response to cold sea- 
surface temperatures at the vortex center (Figure 16). 

The region of maximum upwelling likewise moved outward, 
leaving a broad region of 20C sea-surface temperatures 
between 15 and 45 km and zero mixed-layer depth to a 
radius of 95 km. The dramatic effects of advection were 
also evidenced by significant ocean heat content changes 
(Figure 17). The positive values correspond to heat loss 
regions and the negative values to heat gain regions. 
Vertical and radial velocities produced a redistribution 
of heat by changing the ocean temperature profile (Figure 
18). Because of the ocean's response to the atmospheric 
stress profile, radial currents (u^_) transported heat 
across the model's open boundary at 300 km. Consequently, 
there's a significant averaged heat loss (Table II) in 
the ocean model for this experiment. 

D. ENTRAINMENT AND CONVECTIVE MIXING 

This experiment was similar to the e n tr a inmen t-on ly 

run except heat flux (Q + Q ) was included in (25) and 

E S 

(26). For the initial run, stress partitioning (29) was 

4 

utilized with C = 10 . Radial profiles of mixed-layer 

depth (Figure 19) and temperature (Figure 20) closely 
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Figure 14. Same as Figure 12 except mixed-layer 
depth . 




temperature . 
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Figure 16. Same as Figure 12 except tangential 
component of the wind. 
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Figure 17. Same as Figure 12 except ocean heat 
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Figure 18. Predicted 6- and 12-hour vertical 
temperature profiles at 30 km for 
advection only. 
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Figure 19 . Radial profiles for mixed-layer depth 
at 6 and 18 hours for entrainment and 
convective mixing . 




temperature . 
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resemble the corresponding profiles for entrainment-only 
(Figures 8 and 9). The inclusion of the heat flux term 
produced a slight increase in the amount of deepening and 
cooling compared to entrainment-only; however, the amount 
was less than 2 meters and 0 . 5C at 18 hours. One may 
conclude that entrainment mixing dominates when compared 
to convective mixing for the initial conditions and atmos- 
pheric forcing utilized in the model. 

For this run the oceanic heat content should not be 
conserved because of the latent and sensible heat exchange 
with the atmosphere. Figure (21) depicts radial profiles 
of oceanic heat loss, calculated in the ocean model, and 
heat extraction (H) from the boundar y - lay er model. The 
difference between the profiles represents the loss attri- 
buted to inaccuracies in the numerical solution of (25), 
(26), and (27), and unfortunately is as large as H itself. 
The vertical ocean temperature profiles (Figure 22) were 
similar to those for ent r a i nment -on ly (Figure 11). For 
comparison, Figures (23) and (24) show radial profiles of 
mixed-layer depth and temperature, respectively, for an 
entrainment-convection run in which stress partitioning 

was not used. All available stress (T ) from the boundary- 

a 

layer model was used in the kinetic energy term (G) of 
(25) and (26) to produce changes in mixed-layer depth and 
temperature.. The increased atmospheric forcing was 
evidenced by an increased maximum mixed-layer depth (12 m 
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Figure 21. Predicted 18-hour radial profiles of 
ocean heat loss and accumulated 
from boundary - lay er model. 
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Figure 22. Same as Figure 11 except for entrainment 
and convective mixing. 



51 




Figure 23. Radial profile of mixed-layer depth 
at 18 hours for entrainment and 
convective mixing (no partitioning). 




Figure 24. Same as Figure 23 except mixed-layer 
temperature . 
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deeper) and a lower minimum mixed-layer temperature (.5C 
colder) than for the partitioned run. In both of these 
entrainment-convection runs, as well as the entrainment- 
only run, the amount of sea-surface cooling near the vor- 
tex center was sufficient to cause the radius of maximum 
wind (r_^) to migrate outward. The intensification of the 
atmospheric forcing at larger radii produced a correspond- 
ing increased deepening and cooling of the mixed layer at 
larger radii. Since there was no upwelling in these runs, 
the radial profiles of mixed-layer depth and temperature 
were not modified, and the radius of maximum cooling and 
deepening remained at 20 km. 

In the next set of experiments advection will be in- 
cluded and the effect of variations in the partitioning 
constant (C) will be examined. 

E. ENTRAINMENT, CONVECTION, AND ADVECTION 

This run incorporated all the mechanisms believed to 

be of importance in cooling the ocean beneath a hurricane. 

The total available atmospheric stress (T ) was partitioned 

a 

(29) into that which produced surface turbulence (T^) an< ^ 

that which generated wind-driven currents ( T — T ) . For 

a o 

the initial run the partitioning constant (C) in Equation 

4 

(29) was set equal to 10 . This gave a stress vs. depth 

profile as depicted in Table I. In addition to the kinetic 
energy and flux terms, advective terms were included in the 
solution for mixed-layer depth (26) and temperature (25). 
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The radial advection velocity (u^) was specified according 
to (33) and areas of vertical velocities (w) were calculated 

i 

from (34). Radial profiles of u^ (Figure 25) and w (Figure 
26) were similar in nature to those for advection only 
(Figures 12 and 13). However, since stress partitioning was 
not used for the advec t i on- only experiment, the magnitudes 
of u^ and w were smaller in the current experiment. 

As in previous runs, the hurricane eye wall moved out- 
ward in response to cold sea-surface temperatures near the 
vortex center (Figure 16). The upwelling region (Figure 26) 
likewise moved outward in response to the movement of maxi- 
mum surface stress. The changes in mixed-layer depth and 
temperature produced by entrainment and convective mixing 
(Figures 19 and 20) were significantly mod ified by the out- 
ward migration of the upwelled region. For the entrainment- 
convection mixing case the region of maximum mixed-layer 
deepening was coincident with the region of maximum cooling. 
However, with the addition of advection, the maximum deepen- 
ing (Figure 27) occurred at radii just outside the radius 
of maximum wind (r^), where upwelling was weak and the atmos- 
pheric forcing (G and Q + Q ) was fairly large. At smaller 

F S 

radii the upwelling dominated the effect of mixing, leaving 
very shallow mixed-layer depths. The region of maximum 
ocean cooling (Figure 28) , which remained in the vicinity 
of 15 km, represented an interaction between mixing and ad- 
vection. By decreasing the layer thickness, upwelling 
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Figure 25. 6- and 18-hour radial profiles of ocean 

radial velocity for entrainment, 
convection, and advection (partitioning, 
C=1C) ) . 




Figure 26. Same as Figure 25 except ocean 
vertical velocity. 
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Figure 27. Same as Figure 25 except mixed-layer 
depth . 




Figure 28. Same as Figure 25 except mixed-layer 
temperature . 
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enhanced the effects of atmospheric forcing on the mixed- 
layer cooling (25). Stress partitioning increased the 
percentage of total stress used in the kinetic energy term 
(G) as the mixed-layer depth decreased. There was a corres- 
ponding decrease in the stress used for currents (upwelling) 
which tended to keep the mixed-layer depth from becoming 
zero. Since ultimately it was the value of the constant 
(C) which determined whether mixing or upwelling would have 
the greatest influence on the mixed-layer depth, the entire 
formulation of stress partitioning deserves considerable 
future study . 

Significant ocean heat loss (Figure 29) occurred inside 
30 km where upwelling had the greatest effect on the verti- 
cal temperature profile (Figure 30). At larger radii, where 
entrainment and convection dominated, ocean heat content 
changes reflected the heat extraction (H) by the boundary- 

layer model. The area averaged heat loss (Table II) for 

2 

the solution domain was 2300 cals/cm at 18 hours. As in 

the advection-only case, this loss reflects the transport 

of heat across the model's open boundary at 300 km. 

To test the influence of stress partitioning (29), a 

3 

second experiment was run using a value of C = 5 x 10 . 

With this smaller value, the percentage of total stress 

(T ) used for entrainment mixing was smaller and decreased 
a 

more rapidly with mixed-layer depth and, consequently, the 
percentage of stress used for current generation (T -T^) 
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Figure 29. Predicted 6- and 18-hour radial profiles 
of ocean heat loss for entrainment, 
convection, and advection. 
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Figure 30. Same as Figure 11 except entrainment, 
convection, and advection. 
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t 

was larger and increased more rapidly with mixed-layer 

depth (Table I). Comparing radial profiles of mixed- 

4 

layer depth and temperature for C = 10 [Figures (27) and 

3 

(28)] with those for C = 5 x 10 [Figures (31) and (32)], 

it is apparent that advection by horizontal and vertical 

velocities had a greater influence in the latter case. 

[Note similarity between Figures (31) and (32) and those 

. 3 

for advection only Figures (14) and (15)] . For C = 5 x 10 
upwelling has eliminated the mixed layer in a small region 
near the vortex center and the combined effects of upwell- 
ing and mixing have cooled the sea-surface temperature to 

20C. At larger radii where mixing dominates, there was 

4 

greater deepening and cooling with C = 10 than with 
C = 5 x 10 3 . 

Table II summarizes some significant atmospheric and 

ocean variables. The air-sea interaction was evident in 

all cases where there was pronounced cooling of the ocean 

surface near the vortex center. The lower sea-surface 

temperatures reduced the air-sea temperature difference, 

hence the heat fluxes from the ocean to the boundary layer 

were reduced. The tendency for the model's eye wall to 

seek a region of maximum equivalent potential temperature 

(0 ) caused the radius of maximum wind (r . ) to move outward 

e l 

as the inner surface was cooled (Figure 16) . This outward 

migration maintained a sufficient equivalent potential 

temperature difference (0-0 ) and, consequently, the 

e . e 
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Figure 31. Same as Figure 27 except C=5xl0 . 




3 

Figure 32. Same as Figure 28 except C=5xl0 . 
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maximum wind (v n ) remained fairly constant. Because of 

u . 

1/2 1 

the v r = constant profile, the outward migration of 

u 

r^ produced an unrealistic intensification of the circula- 
tion at large radii. This in turn maintained a constant 
total ocean mass transport with radius which resulted in 
significant amounts of heat being transported out the 
model's 300 km boundary to a region where downwelling would 

presumably have occurred. In order to properly model the 

1/2 

movement of the eye wall the v Q r = constant profile must 

o 

be modified to prevent this unrealistic intensification at 
large radii . 
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TABLE II 



Significant atmospheric and oceanic 
variables for different experiments 
at 18 hours . 





max 

h 

o 

(m) 


min 

T 

s 

(C) 


v e . 

i 

(m/sec ) 


r . 

(km) 


Area Aver- 
age ocean 
heat l°s^ 
(cals/cm ) 


Area 

Average 

H 2 
(cals/cm ) 


Convection 


35 


29 . 5 


60 


21 


252 


240 


Entrainment 


75 


28 . 5 


58 


54 


137 


— 


Advection 
(No Parti- 
tioning) 


32 


20 . 0 


60 


78 


5300 


— — — 


Entrainment 

and 

Convection 














Partition- 

ing 


7 6 ' 


28 . 4 


58 


57 


437 


262 


No parti- 
tioning 


88 


27 . 9 


56 


51 


578 


214 


Entrainment , 
Convection 
and 

Advection 




• 










4 

C = 10 


66 


21.7 


59 


72 


2300 


288 


C = 5 x 10 


3 56 


20 . 0 


59 


66 


1854 


260 
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IV. SUMMARY AND CONCLUSIONS 



The results of this study indicate the feasibility of 
coupling an ocean model to an atmospheric vortex model to 
simulate time-dependent air-sea interaction. Since the 
model was stationary, and the effects of large-scale ocean 
current, turbulent- scale energy dissipation, radiation, 
and internal waves were not included, the magnitude of the 
ocean thermal changes may not be representative. However, 
the relative changes for the various cooling mechanisms 
(convective mixing, entrainment mixing, and advection) were 
apparent. 

First, convective mixing had insignificant effects on 
an ocean thermal structure initialized with 30 meter mixed- 
layer depth, 30C mixed-layer temperature, and below-layer 
gradient of lOC/lOOm. Convective cooling would have been 
more pronounced for shallower layer depths and stronger 
below-layer gradients; however, these were not tested in 
the model. 

Second, entrainment (wind) mixing produced fairly sub- 
stantial mixed-layer deepening and cooling in the region 
of maximum wind stress. The radius of maximum wind res- 
ponded to the colder sea-surface temperatures by migrating 
outward. The effect of entrainment mixing, likewise, 
moved outward producing deeper layer depths and colder sea- 
surface temperatures at larger radii. The depth and 
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temperature changes due to entrainment and convection were 
comparable to those for entrainment only indicating en- 
trainment to be the dominate process. 

Third, the effects of the wind-driven circulation, 
when acting alone, produced the most pronounced changes 
in mixed-layer depth and temperature. • Upwelling inside 
the radius of maximum wind displaced the relatively warm 
mixed-layer water by cold deep-ocean water advected from 
below. The warm sur f a c e - 1 ay er water was advected outward 
by radial currents. Since no significant downwelling 
occurred within the solution domain, heat was lost through 
the model’s outer boundary. 

In reality the various mechanisms studied do not act 
separately but rather interact in some complex fashion to 
produce the thermal changes which observations have shown. 
It appeared to be a question between entrainment mixing 
and advection as to which has the dominate role in cooling 
sea-surface temperatures beneath a hurricane. Advection 
definitely produced significant cooling inside the radius 
of maximum wind but to think of a wind-driven circulation 
with no surface turbulence seems unrealistic. It appeared 
that the interaction of the two processes, with a budget- 
ing of the total atmospheric stress between turbulence 
generation and current generation, was the dominate 
mechanism for sea-surface cooling. 
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